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Abstract. We analyzed observation data collected by the Instrument for the Detection 
of Particles (IDF) on board of DEMETER satellite during the period of total seventeen 
months in 2007 and 2008. In the meantime, the VLF transmitter located at NWC ground 
station was shutdown for seven months and working for total ten months. Our analy- 
sisfor the first time revealed in details the transient properties of the space electrons in- 
duced by the man-made VLF wave emitted by the transmitter at NWC. First, we mapped 
the electron flux distribution and figured out the special range what the NWC belt cov- 
ered. Then we investigated the NWC electron spectrograms in a wide range of Mcllwain 
parameter (up to L^3.0). Finally, we obtained the averaged energy spectrum of the NWC 
electrons within the drift loss-cone, and compared the difference during the observations 
between daytime and nighttime. Our results proved the fact that the VLF emissions from 
NWC transmitter created momentary electron enhancement with fluxes up to 3 orders 
of magnitude. These electrons are distributed in the region of 180° in longitude and 1.6 
~ 1.9 of L shell. In addition, the VLF emission induced either enhancement or loss of 
electrons in higher magnetic shells up to L^3, and the maximum loss was up to 60% 
of the original value. The energy spectra of these electrons revealed that the enhance- 
ment during the NWC daytime are more attenuated than those in the NWC nighttime- 
and that the shape and the cutoff energy of the spectra are also quite different. We will 
present the results of our analysis, compare it with previous studies, and discuss the agree- 
ment of our results with the theory of wave-particle interaction. 



1. Introduction 

More than thirty years ago, it has been learned that the 
VLF wave emitted by transmitters on the ground stations 
can cause precipitation of electrons in the radiation belt, 
and the process obeys the theory of wave-particle interac- 
tion. Subsequently, there have been numerous experimental 
observations and theoretical interpretations on this subject. 
BuUough et aL(1976) studied the ELF/ VLF radiation at the 
altitude of the Ariel sateUite and revealed the phenomenon 
of wave-particle interaction in a geomagnetic conjugate re- 
gion in the southern hemisphere(2 < L < 3). Kimura et 
aL (1983)found strong correlations between the 0.3-6.9keV 
electron fluxes observed by EXOS-B satellite and the 0.3- 
9kHz VLF wave emitted by transmitter at SIPLE ground 
station. 

In the SEEP experiment conducted from May to Septem- 
ber in 1982, by turning the transmitter at the NAA ground 
station on and off at the time scale of seconds, Imhof et 
al (1983) studied the instantaneous correlation between the 
VLF signal emitted from the ground station and the elec- 
tron fluxes observed by the low-orbiting satellite. Inan et 
al. (1985) interpreted the observation of the SEEP exper- 
iment using the theory of particle-wave interaction. They 
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proposed that the precipitation of particles was restricted 
by pitch angle distribution of those particles close to the 
loss cone. 

During the same period, William et al. (1993) studied 
the VLF wave propagation in the D layer of the earth iono- 
sphere, they proposed a multiple-mode three-dimensional 
model. Abel and Thorne (1998) proposed that in the inner 
radiation belt the electron flux reduction mainly induced 
by the VLF wave scattering, which cause the electrons fall 
into the loss cone. Their theoretical calculations indicated 
that the energy of wave-particle resonance decreases with 
higher L values. Furthermore, Richard et al.(2005) studied 
the mechanism of electron acceleration in the outer radiation 
belt. Their studies indicate that these electrons could be ac- 
celerated to MeV energy region by VLF wave with several 
kHz, and the electron fluxes in the observed region could 
rise to 3 orders of magnitude within one to two days. 

The DEMETER sateUite (Detection of Electro-Magnetic 
Emissions Transmitted from Earthquake Regions) launched 
by France in June 2004 is a low-orbiting electro-magnetic 
satellite with detectors onboard to measure space electric 
fields, magnetic fields and high energy particles and so on 
( Sauvaud et al. [2006]; Berthelier et al. [2006]; Parrot et al, 
[2006]). These detectors provide us with the opportunity to 
further study the interactions between the man-made VLF 
wave and particles in the radiation belt. 

By turning on and off the transmitter based on the NPM 
ground station, Inan, Graf et al studied the correlation be- 
tween the ground VLF wave and the electron fiuxes de- 
tected by on-board IDP (Instrument for Detecting Parti- 
cles). The experiments were performed in ten different time 
scales ranging from second to minute, and they found cor- 
relations in 0.1 and 0.2 Hz ON/OFF frequency ( Inan et al. 
[2007]). However, subsequent experiments using the 0.1 Hz 
frequency (i.e. on for 5 sec and off for 5 sec, alternately ) 
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indicated that the correlation rate is only 13.9% ( Graf et 
al. [2009]). Quite recently P. Wang et al.(2010) pointed out, 
the lower correlation rate may be caused by that IDP only 
can measure the electrons with big pitch angles. 

Parrot, Sauvaud, Gamble, et al. studied the relation 
between the VLF emission from the transmitter at NWC 
ground station and the observation by IDP in relatively 
larger regions in space and longer time intervals. There were 
very strong ionosphere disturbance detected near the NWC 
ground station ( Parrot et al. [2007]). Moreover, due to the 
VLF wave emitted by the NWC transmitter, plenty elec- 
trons precipitate into the drift loss cone and drift eastward 
until into the SAA (South Atlantic Anomaly) to covered a 
region with L value between 1.4 1.7 during the local night 
period. The study also reported that the spectrum of local 
electrons has a wing-like structure ( Sauvaud et al. [2008]; 
Gamble et al. [2008]). There are also more relevant stud- 
ies, for example, the numerical simulation for precipitation 
characteristics induced by five ground-based VLF transmit- 
ters ( Kulkarni et al. [2008]), the restriction of VLF wave in 
(18-25 kHz) in driving electron precipitation in inner radia- 
tion belt, as well as the effect of NAA station to electrons 
in outer radiation belt ( Glilverd et al. [2008]; Glilverd et al. 
[2010]). 

Furthermore, recent studied have derived the theoretical 
models of the NWC VLF wave transmission in the iono- 
sphere ( Lehtmen and Inan [2009]), and the different effect 
by ducted or non-ducted VLF wave propagation ( Rodger 
et al. [2010]). A report revealed the increasing of the MF 
component at the height of the DEMETER orbit over the 
VLF transmitters ( Parrot et al. [2009]). It is caused by 
the emissions of the global thunderstorm activity which are 
less attenuated when they cross the ionosphere above these 
locations. 

We analyzed the IDP data during 2007 and 2008, when 
the NWC transmitter was turned off from June 2007 to 
February 2008 and turned on in the rest of the time. NWC 
transmitter has strong radiation power and the very narrow 
bandwidth of the VLF emissionit provided us with the best 
opportunities in studying the mechanism of wave-particle in- 
teractions. In particularbecause the NWC station was kept 
either on or off for a long period, we have enough data to 
make a comprehensive study on the effect of man-made elec- 
tron precipitation at high statistical significant level. 

We analyzed the IDP data during 2007 and 2008, when 
the NWC transmitter was turned off from June 2007 to 
February 2008 and turned on in the rest of the time. NWC 
transmitter has strong radiation power and the very narrow 
bandwidth of the VLF emissionit provided us with the best 
opportunities in studying the mechanism of wave-particle in- 
teractions. In particularbecause the NWC station was kept 
either on or off for a long period, we have enough data to 
make a comprehensive study on the effect of man-made elec- 
tron precipitation at high statistical significant level. 

For the first time, we analyzed the data observed by IDP 
during 7 months of NWC off and 10 months of NWC on. 
By using the (ON - OFF) method, we mapped the instan- 
taneous electron belts, which are induced by the VLF wave 
emitted by NWC transmitter and are termed "NWC elec- 
tron belt", or "NWC belt" in brief. We also quantitatively 
studied the time-averaged characteristics of the NWC elec- 
tron belts, for which we will present the distribution range 
along the longitude and magnetic shell of L parameter, the 
time-averaged energy spectrum of the precipitation elec- 
trons, the variation of the energy spectrum in different L 
shell intervals, and the diurnal variation of the electrons en- 
ergy spectrum. The results of the analysis will be presented 
and discussed. 

2. Statistical analysis of the NWC electron 
flux 

The orbit of the DEMETER satellite is at the height of 
670km (after 2006). It is a quasi-sun-synchronous orbit, 
and the orbital inclination is 98. The satellite is three-axis- 
stabilized. During local daytime the satellite flies downward 



(from north to south) and the descending node is 10:30 at 
noon; during local nighttime it flies upward (from south to 
north) and the ascending node is 22:15 at midnight. The or- 
bital period is 102.86 minutes. The repeated period of orbit 
location is around half of a month. 

The NWC ground station locates in the northwest corner 
of Australia at (21.82 °S, 114.15 °E). Its geomagnetic coor- 
dinate is (31.96 °, 186.4 °).The NWC transmitter works on 
19.8 kHz waveband with very narrow bandwidth of 300Hz, 
and the great emission power of IMW. 

The IDP is placed on flank of the DEMETER satellite 
with field of view (FOV) of 32° , which was pointing towards 
the west during the satellite upward flying at night and to- 
wards the east during the satellite downward flying at day 
time. So in the most case, only those electrons with big pitch 
angle close to 90 degree can be detected by IDP. Figure 1 
shows the location of NWC station and a trajectory of a 
satellite orbit in east side. Figure 2 presents electron spec- 
trogram detected by IDP along with this orbit. It should 
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Figure 1. The location of NWC station and the trajec- 
tory of the orbit No.20520 (upward) 




Figure 2. Electron spectrogram versus UT/LT and L 
shell of orbit No. 20520 (upward) 
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Figure 3. The distribution of average electrons 
flux along logarithm (NL and SL). Energy range: 
206±26.7keV. 



LI ET AL: STUDY OF THE NWC ELECTRONS BELT OBSERVED ON DEMETER SATELLITE 



X - 3 




Figure 4. The electron flux distribution with the energy range 108-411keV for observations in 2008. 
Colors indicate the average flux value within each of (1°,1°) pixel. The simplifled unit of counting 
rate(counts/s) is obtained by multiplying the factor of 1.16(cm^-sr) with the averaged flux value. 



be noted that in Figure 2 there are two wing-like spectrum 
structures at lower energy band around L~1.9 in the north 
and L'^l.G in the south, which correspond to the precipi- 
tated electron belt induced by the NWC VLF wave. These 
will be the focus of this study. 

We mapped the electron flux with the data in the entire 
year 2008. The flux was averaged in each of the l°by 1° pix- 
els along the geographic latitude and longitude. The result 
was shown in Figure 4. The white boxes are the main re- 
gions of the VLF man-made electron belts corresponding to 
the "wing" structure shown in Figure lb). The black lines 
mark the value of the Mcllwain parameter L at the height 
of the DEMETER satellite. Figure 4(a) shows the range of 
the VLF electrons "wing" structure detected during the up- 
ward flying. Whereas during the downward flying, the VLF 
electron belts were detected to move toward east, extending 
all the way to the South Atlantic Anomaly (SAA) region in 
the southern hemisphere, as shown in Figure 4(b). 
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Figure 5. The mean and sigma of double Gaussian fit 
for each months of 2007 and 2008. Circles represent the 
first component and the black squares represents the sec- 
ond component. 



The two white boxed regions in Figure 4 (a) (upward 
data) correspond to the electrons belt near NWC position 
(SL: South Local) and the conjugation position of the NWC 
station (NL: North Local); whereas the three white regions 
in Figure 4 (b) (downward data) correspond to the SE 
(South East) region and NE (North East) region as well 
as the part of SE region extending into SAA region, which 
is belong to SE and marked separately. 

Using upward orbit data, the same mapping have been 
made for each of individual month in 2007 and 2008. We 
find that the electrons precipitation belts in NL and SL re- 
gion vanished from June 2007 to January 2008. So that the 
statistical analysis have been made for each of 24 months 
and both of upwards and downwards data respectively. 

We select the longitude range of A in (110°, 170°) for NL 
and SL regions(170° , 260°)for NE and SE regions, the L 
value in (1.6, 1.8) range, and the electron energy range of 
206±26.7keV for the analysis. The average count rate are 
calculated in each of (AL, AA) with unit size of (0.01, 1°) 
and accumulated to get flux distributions within the selected 
region. 

As examples. Figure 3 displays the result with the data 
in May and June 2007 respectively, which can be fit by 
(a)double Gaussian distribution, and (b) single Gaussian 
distribution in logarithm scale of averaged counting rate. 
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Figure 6. The electron count rate distribution vs. L in 
the area of NL. Longitude: 105° ~ 170°. 
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Table 1. The main area of the NWC belts (wisp range). 





NL 


SL 


NE 


SE 
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1.5~1.9 


1.6~2.1 


1.6~2.2 


1.6~2.1 




1.72 
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110E~180E 
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117E 


125E 
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NWC local time 


16:50~22:30 


16;50~22:30 


21:00~7:40AM 


21;00~7;40AM 



Compared with June, the distribution in May has one 
more component with very wide range of higher count rate, 
and the maximum value can be more than 1000 times of 
the average background value. Figure 5 shows the statis- 
tics result for each of 24 months in the two years. We can 
find that it has two components: a background component 
(meanlzbsigmal) with constant low count rate and a vari- 
able high count rate component (mean2±sigma2). 

This high count rate component was distributed in a wide 
range and more variable. It vanished from June 2007 to Jan- 
uary 2008 which is consistent with the NWC off time, i.e. 
the time when there was no signal of the NWC transmission 
as mentioned in the paper by Zhang et al. (2009). So we can 
confirm that the high count rate component is contributed 
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Figure 7. The electron count rate distribution vs. geo- 
longitude in the area of NL, L: 1.5~1.9. 




by the NWC transmission, because its appearance and dis- 
appearance clearly synchronize with the NWC on and off. 

Based on the ON and OFF status shown in Figure 5, we 
selected "ON data" and "OFF data", respectively. The ON 
data include the observations in ten months from March to 
May in 2007 and March to September in 2008, while the 
OFF data include the observations in seven months from 
July 2007 to January 2008. Since the Solar activity was 
keeping relatively quiet during 2007 and 2008, by subtract- 
ing the background "OFF data" from the "ON data", (ON 
- OFF), we can smooth out short term random fluctuations. 
In addition, we can conduct detailed studies on the basic dy- 
namic characteristics of the NWC electron belts using large 
data sets to get better statistical significance. Since the elec- 
tron fiux varies every month and NWC electrons are moving 
all the time, the obtained results refiect the average statisti- 
cal characteristics. The results will be given in the following 
sections. 



3. Spatial distribution of NWC electrons 

We selected the IDP data with electrons energy in 
91~678keV range, L value of (1.4~3) and geo-longitude 
value covering 180 degree. By (ON - OFF) method, the spa- 
tial distribution of NWC electrons have been investigated 
quantitatively. As an example, results for NL region are 
shown in Figure 6 and Figure 7. The average count rate 
of electrons were calculated for whole given geo-longitude 
range in each value of L~L-|-AL. The results for ON data 
(A), OFF data (B), signal (A-B), and the ratio of signal to 
background R (L) were respectively presented in Figure 6. 
Then The average count rate of electrons were calculated 
for whole given L range in each of geo-longitude interval 
A ~ A -I- AA, the results for different condition including R 
(A) were presented in Figure 7. 

Using the same treatment, the analysis for the other three 
regions of (SL, NE , SE) were also carried out. The results 
of R (L) and R (A) for these regions are presented in Figure 
8 and Figure 9. 



Figure 8. Ratio of signal to background R(L) versus 
L- value in other 3 regions. 
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Figure 9. Ratio of signal to background R() versus geo- 
longitude in other 3 regions. 
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Figure 10. In NL region, the energy spectrum in three 
cases of On, Off, On/off in L range of (a), 1.6 1.65, 
(b), 1.9 1.95, (c), 2.05 2.1; and the ratio of Signal to 
background calculated by R=(On-Off)/Off and present 
in (a2), (b2), (c2). 
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(c)SL (d}SE 

Figure 11. Spectrogram of NWC electrons along with L value in four different regions. ( The data with 
statistical fluctuation only has been ignored in the figures.) 




It is clearly shown that the distributions of NWC elec- 
trons have very complicated structures along the longitude 
and the magnetic shell. In SE regionthe effect from NWC 
VLF wave expands up to 300° , which is deeply into S AA re- 
gion due to the eastward drift of NWC electrons. With the 
high background electron flux in SAA, the NWC electrons 
still have S/B by factor of more than f~2 times. The ef- 
fect of expanding to the radiation belt along with magnetic 
shell up to L~3 can also not to be ignored. These suffi- 



ciently indicate that the VLF emission of NWC transmitter 
powerfully influence the flux of the space electrons. 

We deflne the range of drift lose cone for NWC electrons 
belts according to R (L) and R (A) distribution, and give 
the name as "NWC electron belt", or "NWC belt" in brief, 
which are listed in Table f . 

It's worth mentioning that the data include both upward 
data (night for the satellite) and downward data (day for the 
satellite), so overlapping range exists along the longitude di- 
rection. The longitude ranges in Table 1 correspond to the 
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local nighttime of the NWC station from 16:50 to 7:40am. 
Since the increasing of the NWC electron flux can only be 
obviously observed during the NWC night, the belts cover 
half the Earth with 180 degree along the longitude. The 
related ionosphere effect will be discussed in the subsequent 
section 4.2. 

4. Energy spectrum of NWC electrons 

We investigated the energy spectrum of NWC electrons 
and its variation in various situations. The spectra were 
accumulated by On, Off, and (On-Off) respectively. The 
selected region is basically according to Table 1. 

4.1. Variation of energy spectrum along with L shell 
value 

According to the results shown in Section 3, we noted 
that, at a certain altitude, the flux of NWC electrons var- 
ied with L value of the satellite passed. The same effect 
exists for energy spectrum. As an example, we presented 
the result for NL region in Figure 9. The differential energy 
spectrum in three narrow range of L shell are presented for 
ON, OFF, and (ON - OFF) status respectively. We firstly 
revealed that: (a). Within NWC belt (L: 1.6 - 1.65), flux of 
NWC electrons has a peak centered at 220keV with S/B of 
35 times; (b) For higher L region of 1.9 ~ 1.95, the electron 
enhancement reduced and dominated in 300keV IMeV, 
which reach to a maximum value of S/B about 1.3 times 
near 600keV; (c) For region in inner radiation belt L:2.05 
~ 2.1, a slot region is formed by the NWC electrons be- 
low 700keV, with S/B value reaching -60% near the peak of 
300keV. In all three situations, no NWC particles exist with 
energy above IMeV except the statistical fluctuation. 

In order to accurately investigate the evolvement of NWC 
electrons energy spectrum along with L shell value, we per- 
formed further study of spectrogram for the L shell with 
value from 1.3 to 3 in step size of 0.05 for the four NWC 
electrons belts, which have the same longitude range with 
table 1. The obtained results of enhanced electron fluxes 
from ON-OFF value are shown in Figure 11 and that of lost 
electrons shown in Figure 12. 

From Figure 11 and 12, we can find the following charac- 
teristics of NWC electrons energy spectrum: 

1) , In all 4 regions with electrons precipitation, the spec- 
tra show obvious "wing" structures in relative lower L shell 
range. The main part of "wing" structures correspond to 
NWC electrons beltswhich have been defined in Table 1. 

2) , In Local regions, there are obvious slot effect caused by 
electron loss in the energy range of 100 ~ 200keV. The elec- 
tron loss is most obvious in north region (NL) corresponding 
to L shell value in 1.95 ~ 2.5. The spectrum shows irregular 
structure in south region (SL) corresponding to inner radi- 
ation belt with L i 2.2. 

3In Eastward regions, there are no obvious slot effect of 
NWC electrons. The enhancement in NE and SE region 
with energy j300keV extend to L ~ 2.5 and even to L ~ 2.9 
respectively. The fluxes of lost electrons are rather lower, the 
distributions are mainly in 300 ~ 400keV range and slightly 
different between in SE and in NE. 

Related to the physical mechanism, our understanding is 
that, the "wing" structures should be most likely caused by 
local interactions between original 19.8 kHz of VLF wave 
and space electrons, which induce the electrons pitch an- 
gle scattering and precipitation into loss cones. The origi- 
nal VLF wave, most probably, come from NWC transmitter 
by direct ducted propagation. And more electrons covered 
wider L shell range in Eastward are the results from drift of 
precipitation electrons along with longitude and magnetic 
shell, as well as from more electron resource provided by 
SAA and inner radiation belt. 

There are two more factors need to be mentioned here. 
First, our analysis uses a geographic longitude coordinates 



rather than the geomagnetic longitude coordinates. It may 
cause the difference between north and south, east and west, 
for example, the covered L range of "wing" feature are more 
similar between NL and SE, as well between SL and NE. 
Another factor is the "east-west" effect of electron fluxes 
caused by opposite pointing of IDP aperture during upward 
and downward of the orbit. The difference is usually equal 
to around 16% in mid-latitudes of northern hemisphere NL 
(0 ~ 28'="''=) (Li XQ et al., 2010). The effect for NWC belts 
in higher latitude can be roughly estimated from Figure 5. 
The differences of background fluxes between two regions 
could reach to two times or more. Similar effect for NWC 
electrons can also be seen in the figure. 

4.2. Differential energy spectrum of NWC electrons 
and "Day-Night" effect 

By (On-Off) method, we obtained total differential en- 
ergy spectra of NWC electrons, which are accumulated from 
each of 4 regions defined in Table 1. Figure 13(a) presents 
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Figure 13. Differential energy spectrum S(E) of NWC 
electrons in drift loss cone, (a): During NWC night time, 
(b):During NWC day time. 
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Figure 14. The relation between the electron flux and 
the angular distance to NWC. 
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Table 2. Characteristics of energy spectrum of NWC electron belts during NWC nighttime, the values in parenthesis are errors. 



Area 


Average Energy 


Integrate influence 


Integrated Flux 




(keV) 


er g 1 cm?" 1 s 1 sr) 


(electrons /cm^/s / sr) 


NL 


167.0 (11.1) 


3.59 (0.12) 


1344.1 (43.7) 


SL 


169.8 (12.8) 


2.68 (0.10) 


986.8 (36.8) 


NE 


176.1 (7.4) 


7.60 (0.16) 


2696.2 (54.4) 


SE 


185.2 (8.5) 


7.54 (0.17) 


2543.5 (58.0) 



Table 3. Characteristics of energy spectrum of NWC electron belts (during NWC daytime), the values in parenthesis are errors. 



Area 


Average Energy 
(keV) 


Integrate influence 
(lO^^erg/crn^/s/sr) 


Integrated Flux 
(electrons / cm? / s/sr) 


Flux ratio 
(night/day) 


NL 


157.7 (39.7) 


0.096 (0.012) 


38.0 (4.8) 


35.4 


SL 


161.0 (31.7) 


0.601 (0.065) 


233.7 (20.8) 


4.2 


NE 


378.6 


0.071 (0.032) 


11.7 (5.9) 


231.1 


SE 


203.8 (109.5) 


0.063 (0.016) 


19.4 (5.4) 


130.9 



the spectra related to NWC local night time using upward 
data for SL and NL region, and downward data for SE and 
NE region. Figure 13(b) is the spectra related to NWC lo- 
cal daytime using downward data for SL and NL region, 
and upward data for SE and NE region. For each spectrum, 
the average energy, the integrated energy density and the 
integrated flux density, as well as the differences of these 
physical parameters between NWC day time and night time 
have been derived and listed in Table 2 and Table 3, respec- 
tively. 

Here we give brief discussion about the result. 

1) , During NWC local night time, spectrum feature are dif- 
ferent between local and eastward region. The spectrum in 
SL and NL region has almost the same feature, with the 
same peak below 200keV and the same cutoff energy 5 at ~ 
500keV. There is enhancement of electrons flux in Eastward 
region and the cutoff energy increase to 800keV. This phe- 
nomenon indicates the existence of continuous process of the 
electron precipitation and being accelerated during drift. 

2) , During NWC Local daytime, we still can see the NWC 
electrons, but the fluxes of electrons are rather weak, and the 
spectrum features arc quite different with each other. The 
cutoff energy of NL and SL region are increased to 600keV 
and 700keV respectively. The energy spectrum takes on 
anti-symmetric structure, the flux in SL region is strongest 
and NL is weakest. The spectra in eastward region show 
more absorption of electrons on their way of drift, so that 
only the Electrons with E > 250keV in NE region can be ob- 
served. As we know that the average energy of wave-particle 
resonance is ~ 170 koV for 19.8 kHz VLF wave, so it at least 
implicates that the original 19.8 kHz of VLF wave almost 
can not reach to NE region due to ionosphere absorption 
during daytime. 

3) , We find out that the integrated electrons fluxes in day- 
time observation can be well fitted by a quadratic function 
versus the angle distance of D, which are calculated from 
coordinate of central position for each belt related to the 

coordinate of NWC ground station. The result is presented 
in figure 14, and formula Iwhere D denotes the angle dis- 
tance with unit of radian. 



55 68 2 
F = — ^2 2.37electrons/cm /s/sr (1) 

It provided the evidence that the wisp part of NWC elec- 
trons are induced by original 19.8 kHz of VLF emission 
which are directly propagated from NWC ground station. 
It play major contributions to NWC electrons during NWC 
night time, and the strong ionosphere absorption during 
NWC day time. This is consistent with what we discussed 
in the previous subsection. 
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Figure 15. The pitch angle distribution of NWC elec- 
trons in NL and SL region during NWC night time. 



5. Verification of wave-particle theory from 
the pitch angle distribution 

IDF has no ability to identify the pitch angle of incident 
particles. So, the direction between the axis of IDF aperture 
and the local magnetic field is taken as the pitch angle and 
recorded into the database. For each of four NWC electron 
belts, we make such kind of pitch angle distribution, as ex- 
ample, the results for NWC night time in NL and SL region 
are presented in Figure 14. The actual distribution has to be 
expanded by a normalized Gaussian function with 32° FOV, 
here we take 16° as 3cr, to obtain the pitch angle response 




^^3.5 240 24,5 250 

Equatorial Pitcli Angle a^[degi 

Figure 16. Local pitch angle diffusion coefficient as a 
function of equatorial pitch angle for given L and E. The 
area between two vertical lines corresponds to electrons 
with pitch angle ranging from 23.5° to 25° at equator 
which is 70° to 90° at L=1.8 observed by DEMETER. 
Electron's kinetic energy E=220KeV. 
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function of the observational incident electrons. Analysis 
shows that the distributions are quite different in different 
belt regions, as well as between NWC day and NWC night. 
But in all of the results, measurable electron pitch angle are 
distributes mostly within 60° ~ ffO°. 

In the inner radiation belt, based on resonant interac- 
tion after ignoring the energy diffusionthe quasi-linear dif- 
fusion equation (Fokker-Planck equation) can be written 
as{Melrose et al. [1980]): 

-E7 = -■ — {Dec stna ) (2 
ot sma oa oa 

Where f(Q,E,L) is density function in phase space, which 
depends on local pitch angle of a, kinetic energy of E and Mc 
Ilwain parameter of L. Daa is local pitch angle diffusion co- 
efficient, which contains all information about wave-particle 
interaction. As for NWC electrons with a given value of 
L shell and kinetic energy of E, the value of local diffusion 
coefficient Daa rely on the pitch angle distribution. 

Local diffusion coefficient of resonant interaction between 
electrons and VLF wave can be written as Summers et al. 
[2005]: 

1 A l2 n — £i££££\2 I dxj i 

J _ IT He 1 AO \^ y.^ ^ ^dyi\ r . Xi — Xm 

- 2T{E + W^^ 5x\Pcosa-^\ '"^^^"^^^ 

(3) 

Where /3= v/c, B is Earth's magnetic field; Xm and are 
reduced parameter defined as Xm —=(jJm/^e, Sx =Soj/Qe; Xi 
and x/i are reduced variable with Xi =uJi/Q,e, Hi =cki/Qe, re- 
spectively. Assuming pitch angle diffusion is occurred only 
within equatorial plane and is valid for purely field-aligned 
wave propagation, we obtain equatorial pitch angle diffusion 
coefficient. 

According to the DEMETER observation, to calculate 
equatorial pitch angle diffusion coefficient, parameters are 
chosen as follows: uj,n = 19.8kHz as frequency of maximum 
wave power, Slo= 150Hz and wave amplitude Ab =10pT. 
Equatorial magnetic field is given by dipolc model B = 
3.11* IOVL^T and equatorial plasma density No = 880* 
(2/L)'*cm^'* which is within the range given by ( Angerami 
et al. [1964]; Inan et al. [1984]). 

From Figure 16, we can clearly find that electrons pitch 
angle is just within the FOV range or nearby of IDP, where 
electrons participate the interaction with VLF wave trans- 
mitted by NWC. As a result these electrons are observed by 
DEMETER after diffusion process. The reason for the ap- 
pearance of "NWC man-made radiation belt" thus is prob- 
ably explained by pitch angle diffusion based on resonant 
interaction. 

6. Summary and discussion 

Based on DEMETER IDP observations at the altitude of 
670 km, the characteristics of the electrons precipitation, in- 
duced by VLF wave emission from NWC station, have been 
firstly studied in detail and quantitatively. With the long 
term observation data related to NWC on and off, the back- 
ground can be eliminated by (On-Off) method with high 
statistics. For the first time, the space distributions, the 
spectrogram via L shells, and the differential energy spectra 
of NWC electrons have been obtained. 

The results based on NWC night data indicate that, at 
670 km height, NWC VLF wave induced huge enhancement 
of electrons flux in broad space range including the belts 
with "wing" structure. But it also causes the electrons fiux 
reduction which happened in lower energy band and at some 



higher L regions. The electrons precipitation, whatever en- 
hancement or reduction, all can be explained by pitch angle 
scattering of the electrons, which caused by wave-particle 
resonance interactions, and possibly by the electrons drift 
along the longitude. This phenomenon is a kind of dynamic 
balanced process of the electrons flux variations. It depends 
on which process is dominated for the electrons either en- 
hanced locally from higher L shell or moving down to the 
lower L shell. 

Used the data from both of NWC daytime and night time, 
respectively, for the first time, the energy spectra are ana- 
lyzed for each of 4 belt regions. It has been shown that 
there are still electrons enhancements during NWC day- 
time. But the fluxes are rather weak comparing the results 
in NWC nighttime, even the electrons in lower energy band 
are all disappear in eastward region. It's already well known 
that this is caused by VLF wave absorption of ionosphere 
at daytime. Interestingly, this effect presents more regu- 
larity between the electron flux and the distance of VLF 
wave propagations, which shown inverse-quadratic function 
relation (Fig. 14 and equation 1). This easily links with 
the fact that the VLF energy density attenuates by distance 
squared and then requires the linearity relation between the 
flux of the precipitated electrons and the energy density of 
the VLF wave. This analysis provides reference for studying 
2the ionosphere transmission characteristics of VLF wave. 
■I We may study more in further work. 
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